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D. Abstract. We find the limit distributions for a spectrum of a sys- 

f"*^ ' tem of n particles governed by a fc-body interaction. The hamil- 

tonian of this system is modelled by a Gaussian random matrix. 

We show that the limit distribution is a g-deformed Gaussian 
distribution with the deformation parameter q depending on the 
' value of the fraction . The family of q-deformed Gaussian dis- 

■ tributions include the Gaussian distribution and the semicircular 

\ law; therefore our result is a generalization of the results of Wigner 

|Wig67[|Wig57| , Mon and French [MF75j . 
^t" ' PACS: 05.40.+j 

O ! 

O ■ 1. Introduction 

1.1. The A;— body interactions. We say that the interaction of n par- 
\ tides is a &-body interaction if it can be treated as a sum of indepen- 
dent interactions, each taking place in a group of k particles (the groups 
of course need not to be disjoint). The integer number k is called a 

>- . rank of the interaction. 

^ \ From basic principles of physics we would expect that the fundamen- 

tal interactions should take place only in pairs of particles; by taking 
into account higher order interactions between carriers of interaction 
we see that also k-bodj interactions are possible if A; is a small integer. 

Perhaps even better examples of a fc-body interaction we obtain by 
considering models of complex physical systems which were simplified 
by disregarding some components. In such models effective hamiltoni- 
ans contain also A;-body interactions (A; > 2) in order to preserve the 
effects connected to the disregarded components. 

Nevertheless we would expect that k, the rank of the interaction 
should be relatively small. 

1.2. Random matrix models for complex quantum systems. 

1.2.1. Wigner 's model. The first and the simplest random matrix model 
for complex quantum systems such as atomic nuclei was proposed by 
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Wigner |Wig67| . In this model a hamiltonian of the system is repre- 
sented by a large hermitian matrix (<hj)i<i,j<N which entries are com- 
plex Gaussian random variables. This model can be heuristically jus- 
tified as follows: in an sufficiently complex physical system the matrix 
elements of the hamiltonian should be very complicated as well and 
therefore can be regarded as random. 

In the Wigner's model the matrix element between any two physical 
states can take nonzero values, what for quantum systems consisting of 
n parts means that the system is governed by an n-body interaction, 
what in the light of the preceding discussion is physically doubtful. 

1.2.2. Random matrix model for k -body interaction. The following more 
realistic model was proposed [BF71al IBF71bl IFW701 IFW711 IWF72j in 
which a &-body interaction hamiltonian of n particles is modeled by a 
random matrix. 

We consider a quantum system of n distinguishable particles; the 
Hilbert space H of the system is a tensor product of the Hilbert spaces 
Hi (1 < i < n) corresponding to particles. In fact as "particles" we 
can take also quantum statistical objects such as collective excitations, 
holes, etc. 

We assume that the rank of the effective interaction in our system is 
equal to k and therefore the hamiltonian H of the system is a sum of 
hamiltonians Ha of the /c-particle subsystems, 

A 

The sum is taken over all sets A C {1, . . . , n} which have exactly k 
elements. 

Due to the factorization of the Hilbert space H = Ha ® Ha' where 
Ha = ^S>ieA^-^ ft A' = ^S>i^A^-i we can wr he each hamiltonian Ha as 
Ha = H A (g) where H A : Ha — ► Ha is a self adjoint operator and 
1^ , : Ha 1 — > Ha 1 is the identity. Suppose that each of the particles 

k 

has s possible states and therefore Hi = C s and Ha = C s ; we see that 
H A can be viewed as a hermitian matrix with s k rows and columns. 

Similarly as in the Wigner's model we shall assume that H A = (%,•) 
is a hermitian random matrix, i.e. afj (1 < i < j < s k ) are com- 
plex Gaussian random variables with the mean and the covariance 
^( a tj a k.i) = E( a tj a tk) = fnjj; 8ii5jk- We assume that the entries of 
different hamiltonians Ha are independent. 

The above considered random matrices are related to Gaussian uni- 
tary ensemble, but it is easy to write a version which is related to 
Gaussian orthogonal or symplectic ensemble. 
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1.3. Overview of the article. The goal of this paper is to investigate 
the limit distribution of the above spectra when n tends to infity and 
the rank of the interaction k(n) depends in a certain way on n. 

A classical result of Wigner |Wig67[ |Wig57| states that if k{n) = n 
then the distribution of the eigenvalues converges to the semicircular 
distribution (see Fig. [1]). On the other hand Mon and French [MF75] 
showed if k is constant and n converges to infinity then the distribution 
of eigenvalues converges to the Gaussian distribution (see Fig. E]). 

In this article we show that in the intermediate cases when 1 <C 
k(n) <C n the limit distributions are given by so-called g-deformed 
Gaussian distributions. 



2. The g— deformed Gaussian random variables 

For overview articles on g-deformed commutation relations and g- 
deforme d Gaussian vari ables we refer to [FB701 IBS9T1 IBS941 IBKS971 



lvLM95t ISniOlaj ISni01b[ Spe92|. In this article we will consider only 



one g-deformed Gaussian variable, what simplifies the discussion sig- 
nificantly. 

Let us consider a real number — 1 < q < 1. One says that an 
operator a acting on a some Hilbert space K, and its adjoint a* fulfil 
the g-deformed commutation relations if 

aa* — qa*a = 1, 

where 1 is the identity operator. 

Suppose that there exists a unital vector Q 6 /C such that 

aVl = 0. 

A vector with such a property is called a vacuum. In such a setup 
we can introduced a (non-commutative) expectation value t(X) = 
(Q,XQ). 

By the g-deformed Gaussian random variable we mean the operator 
a + a*. 

The g-deformed Gaussian distribution is a compactly supported 
probabilistic measure v q on the real line with a property that for each 
natural number n we have 

x n dv q (x) = r[(a + a*) n ] = (a + a*) n O). 

It can be proven |Sze26j that for — 1 < g < 1 this measure is supported 
on the interval — ?=, —3= and the density of this measure is given 
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Figure 1. Wigner semicircular law, corresponding to 
q = 0. 

by 

1 oo 

uJdx) = - JT^sm9T](l -q n )\l -q n e 2ie \ 2 dx, 

n=l 

where 9 G [0, 7r] is defined by x = cos 6 1 . 

It turns out that for q = 1 the g-deformed Gaussian distribution 
coincides with the standard Gaussian distribution (see Fig. |6]). For 
q = the g-deformed Gaussian distribution coincides with the semi- 
circular distribution of Wigner (see Fig. [p. The g-deformed Gaussian 
distributions for intermediate values of the deformation paramater q 
are presented on the Figures [2] [51 



3. The distribution of eigenvalues 



The connection of the random matrix models considered in Intro- 
duction with g-deformed Gaussia n variables was found by the author 



in the following theorem [SniOlb 



fc(ra) 

distribution of the eigenvalues of the hamiltonians H is the q-deformed 



= c, where < c < oo then the limit 



Theorem 1. If lim r 

'the eig 

Gaussian distribution with the parameter q given by 



q = exp 



1 

s z 



where s is the number of one-particle states. 
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Figure 2. The g-deformed distribution for q = 0.2. 
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FIGURE 3. The g-deformed distribution for q = 0.4. 



It should not be a surprise to see that the above theorem contains 
results of Wigner (in the model considered by him we have k = n and 
lim^oo -^U = oo, q = exp[— oo] = 0) and of Mon and French (k is 
constant, therefore c = and hence q — 1). 

We see that the square root from the number of particles is the scale 
of the rank of the interaction in which the passage from the Gauss law 
to the semicircular law occurs. 
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Figure 4. The g-deformed distribution for q = 0.6. 
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Figure 5. The g-deformed distribution for q = 0.8. 
4. Final remarks 



Particles considered in this article were distinguishable. The next 
step would be to replace them by fermion or bosons. In the fermionic 
case it would mean for example that s, the dimension of a one-particle 
Hilbert space must be a function of n, the number of particles. Simple 
arguments show that if s(n) tends to infinity faster than linearly, then 
the most of the states are not occupied and Pauli exclusion principle 
does not affect our calculations. In this case we can take s = oo and 
therefore 
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Figure 6. The Gaussian distribution, corresponding to 
q=l. 
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